Spontaneous currents and half-integer Shapiro steps in 
Superconductor-Ferromagnet-Superconductor 0-n Josephson junctions 
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We study Superconductor-Ferromagnet-Superconductor (Nb-Cuo.47Nio.53-Nb) Josephson junctions 
with spatial variations in the barrier thickness. Critical current vs. magnetic flux diffraction pat- 
terns indicate that the critical current density changes sign along the width of the junctions, creat- 
ing interfaces between and 7r junction regions around which spontaneous currents can circulate. 
Shapiro steps observed at half-integer Josephson voltages can be explained by the phase-locking of 
the spontaneous circulating currents to the applied rf modulation. 



In the past few years considerable attention has 
been directed toward the understanding and realiza- 
tion of 7r Josephson junctions lj. Transitions be- 
tween and 7T junction states were demonstrated 
in Superconductor-Ferromagnet-Superconductor (SFS) 
junctions as a function of temperature |2j and barrier 
thickness |3j, and in mesoscopic Superconductor-Normal 
metal-Superconductor (SNS) junctions as a function of 
the barrier thermalization voltage Q. In both systems, 
the first-order Josephson supercurrent vanishes at the 0- 
7r transition, making the transition region feasible for the 
observation of second-order Josephson tunneling charac- 
terized by a sin(2<^>) component in the current-phase re- 
lation (CPR). Period doubling in the critical current I c 
vs. magnetic flux $ patterns of dc SQUIDs incorporat- 
ing two mesoscopic SNS junctions, one biased at the 0-tt 
transition, has been attributed to a second-order Joseph- 
son component 0. In SFS junctions, Shapiro steps at 
half-integer Josephson voltages were reported close to a 
temperature minimum in I c , consistent with a sin(20) 
Josephson component [||. However, direct measurements 
of the CPR in uniform SFS junctions revealed no dis- 
cernible sin(2(/>) term jjj. 

In this paper, we present experimental evidence that 
half-integer Shapiro steps in the current-voltage charac- 
teristics can occur in SFS junctions with a non-uniform 
critical current density without the need for an intrinsic 
sin(2</>) component. In our junctions, we observe half- 
integer Shapiro steps, but only in a narrow temperature 
range close to the temperature at which a deep but finite 
minimum occurs in the zero- field critical current. In this 
regime, we find that part of the junction is in the state 
and part is in the tt state with the net critical currents 
of the two regions comparable in magnitude, and that 
the energy of the junction is minimized by generation of 
a spontaneous circulating current. This current can cou- 
ple to the applied microwave field, producing half-integer 
Shapiro steps. A similar effect has been observed in an 
analogous system, a nearly-symmetric dc SQUID with an 



applied magnetic flux of i<3?o ($o = h/2e) jgj. 

A 7r Josephson junction is characterized by a nega- 
tive critical current. The mechanism of the it state in 
SFS Josephson junctions is the spatial oscillation of the 
proximity-induced order parameter inside the ferromag- 
netic barrier which arises from its exchange field [3] . The 
critical current density J c is predicted to oscillate and 
decay with the barrier thickness d according to fl?T | 
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where £pi is the decay length and 2tt^f2 is the oscillation 
period of the order parameter. This expression is valid 
for d ^> The lengths £pi and £>2 can be extracted 
by fitting the measured J c {d) to Eq. Q). In our junctions 
at T = 4.2 K, £pi « 1.3 nm and ~ 3.7 nm, and the 
first two nodes of J c occur for d 11 nm and d 22 
nm. The lengths £fi and £f2 vary with temperature 
according to 
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where E ex is the ferromagnet exchange energy and 
£fi,F2(0) are the values at zero temperature. By using a 
weakly-ferromagnetic alloy Cuo.47Nio.53 with Tcurie ~ 60 
K as a barrier material, temperature changes in the 1-4 
K range have a significant effect on the suppression and 
modulation of the induced pair correlations in the ferro- 
magnetic interlayer, allowing us to tune through the 0-tt 
transition by changing the temperature. 

The SFS junctions were patterned by optical lithog- 
raphy. Base and counterelectrode superconducting lay- 
ers were dc-sputtered Nb with thicknesses 100 nm and 
240 nm respectively, separated by an 11 nm barrier layer 
of rf-sputtered Cuo.47Nio.53 and a 20-30 nm layer of dc- 
sputtered Cu. The ferromagnetic layer thickness is cho- 
sen near the first 0-7T transition thickness, while the Cu 
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layer protects the barrier during processing. Junction 
sizes were 4 /zm x 4 /im or 10 /im x 10 /im, denned by 
a window in an insulating SiO layer deposited on top of 
the CuNi/Cu barrier. Because SFS junctions have small 
I c Rn products ~ 1-100 nV, a commercial dc SQUID with 
a standard resistor R st 10 mfi is used as a potentiome- 
ter with sensitivity 1 pV to perform transport measure- 
ments. This translates into a critical current measure- 
ment resolution of 100 nA. A uniform magnetic field up 
to 100 G can be applied through the junction barrier from 
a solenoid coil to measure the magnetic field dependence 
of the critical current, and an rf current component at 
frequencies f r f = 0.3-1.3 MHz is superimposed with the 
dc bias current to observe ac-induced Shapiro steps in 
the current- voltage characteristics. 

Measurements of the critical current I c vs. the ap- 
plied magnetic flux $ threading the junction barrier re- 
veal that the critical current distribution is often not uni- 
form across the junction. Figure db) shows a series of 
Ic{&) curves in the temperature range 1.4-4.2 K for a 
10 /im x 10 /im junction. At T — 4.2 K, I c ($) has a 
Fraunhofer-like shape but with non-vanishing supercur- 
rents at the side minima. This can occur in a junction 
with a localized region of high critical current density. 
In the temperature interval 1.4-1.9 K a minimum in the 
critical current is observed at zero field, indicating that 
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FIG. 2: (a) Calculated temperature variation of the critical 
current for junction sections 1 and 2 and for the entire junc- 
tion. The thin region remains in the state at all tempera- 
tures while the thick region crosses from the state into the n 
state, (b) Measured critical current vs. temperature at zero 
applied magnetic field dips but remains finite. 
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FIG. 1: (a) Stepped ferromagnetic barrier deduced from crit- 
ical current vs. applied magnetic flux measurements, (b) 
Diffraction patterns at a series of temperatures showing de- 
viations from Fraunhofer behavior at low temperatures, (c) 
Simulated diffraction patterns using the deduced ferromagnet 
barrier profile. 



regions of opposite polarity critical current density exist 
in the junction. 

The temperature evolution of the / c ($) patterns indi- 
cates that some fraction of the junction width makes a 
transition from the state to the ir state as the temper- 
ature is lowered, while the remaining part stays in the 
state. The critical current non-uniformity likely arises 
from spatial variations in the barrier thickness across 
the junction but could also be caused by inhomogeneities 
in the ferromagnet exchange energy or by variations in 
the S-F interface transparency. Figure d a ) shows the 
step barrier geometry deduced by fitting the measured 
diffraction patterns in Fig. db) in the short junction 
approximation in which magnetic fields from the tunnel- 
ing current are neglected. We obtain good agreement as 
demonstrated in Fig. dc). It is not surprising that the 
small 6 A step has such a dramatic effect on the diffrac- 
tion patterns since close to the 0-7T transition the critical 
current density in our junctions changes by 1000 A/cm 2 
per 1 nm change in the barrier thickness. For the barrier 
profile in Fig. da) and the experimental data for J c {d), 
we use Eqs. d an d © to calculate the temperature 
dependences of the zero-field critical currents I c \ of the 
thin narrow region, I C 2 of the wide thick region, and I c , 
the total junction critical current. These are plotted in 
Fig. d a )- We see that I c \ is relatively constant while 
I C 2 decreases and changes sign at T w 2.1 K, causing the 
I c to vanish at T n Q w 1.55 K. However, measurements 
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plotted in Fig. |2h) show that I C (Q — 0) does not go 
fully to zero at T v q, instead reaching a minimum value 
of « 10 jit A. 

Shapiro steps induced in zero applied magnetic field at 
drive frequency f r f are also anomalous, exhibiting not 
only the usual steps at integer multiples of the Joseph- 
son voltage V n = n(hf r f /2e), but also steps at half- 
integer Josephson voltages such as V1/2 = hf r f/4e and 
1/3/2 = 3/i/ r //4e when the temperature is close to T v o- 
An example is given in Fig. Efa) for f r f = 1.3 MHz. Fig- 
urc[3f b) shows the maximum value of the (V=0) Joseph- 
son supercurrent, and the maximum amplitudes of the 
n=| and the n=l Shapiro steps obtained by adjusting 
the rf-amplitude (the n th step amplitudes vary with rf 
voltage V r f according to the corresponding Bessel func- 
tions J n (2eV r f /hf r f) as expected). At temperatures far 
from T„.o, only integer Shapiro steps were observed. Half- 
integer steps appear near the minimum in the critical cur- 
rent in a temperature range of width ~ 60 mK. Note that 
the minimum in the critical current measured in Fig. EJb) 
is about 25 mK lower than that in Fig. EJb) as a result of 
room-temperature annealing of the ferromagnetic barrier 
duringthe 3 days between when the measurements were 
made |2(. In other junctions in which the critical current 
does vanish at T^o, half-integer Shapiro steps were not 
observed at any temperatures. 

The non-vanishing critical current at T^o and half- 
integer Shapiro steps could result from a sin(2</>) com- 
ponent present in the current-phase relation near the 0-n 
transition. However, we believe that they are more read- 
ily explained by self-field effects that must be taken into 
account in finite width 0-7r junctions. At temperatures 
close to T„o, spontaneous currents circulate around inter- 
faces between and 7r regions to lower the total energy 
of the system 0, 0. These circulating currents gen- 
erate magnetic flux through the junction that prevent 
the total critical current from vanishing at any applied 
magnetic field. In addition, they can resonate with an 
rf bias current applied at twice the Josephson frequency, 
resulting in half-integer Shapiro steps. A detailed calcu- 
lation of the spatial distribution of spontaneous currents 
and the net critical current in this regime requires a self- 
consistent solution of the Sine-Gordon equations in the 
junction. Computation of the amplitudes of the half- 
integer Shapiro steps further requires numerical simula- 
tions of the junction phase dynamics. However, we can 
understand the onset and effects of the spontaneous cur- 
rents in a 0-7T junction in analogy to a dc SQUID with 
finite geometric inductance. 

We consider a dc SQUID with a junction of criti- 
cal current I c q > and a 7r junction of critical current 
I C7r < in a loop of inductance L. In such a SQUID, 
phase coherence around the SQUID loop precludes both 
junctions being in their low energy states and the junc- 
tion phases 4>o and <p„ depend on the inductance param- 
eter (3 L = 2irLI c /$ , where I c = (\I c0 \ + \I C7r \)/2, and 



the critical current asymmetry a = (|/ c o| — \I c -n\)/2I c . 
For /3l < 2a/(l — a 2 ), it is energetically favorable to 
switch the phase of the junction with the smaller criti- 
cal current magnitude into its high energy state, e. g. 
if |7 c o| > \Icw\i then <f>Q = 4>„ = 0. Under the influ- 
ence of a harmonic ac-drive, </>o and (j)^ wind in phase 
and only integer Shapiro steps can be observed. When 
Pl > 2a/ (1 — a 2 ), the SQUID energy is lowered by gener- 
ation of a spontaneous circulating current and 4>o 7^ <Ptt- 
In this regime, the phase differences (f>o and cjj^ are no 
longer synchronized. Hence, the spontaneous circulating 
current J = (I c / (3l){4>o — 4>tt) can phase-lock to the driv- 
ing frequency in such a way that it switches direction 
an even or odd number of times during each period of 
the drive signal, corresponding to integer or half-integer 
Shapiro steps. This phenomenon has previously been 
studied in an equivalent system, an ordinary dc SQUID 
with an applied flux of j^o- Measurements and simula- 
tions showed half-integer Shapiro steps with amplitudes 
that increase as the SQUID inductance is increased and 
the critical current asymmetry is decreased 8J. 

In a 0-7r junction, the existence of spontaneous cir- 
culating currents, and hence half-integer Shapiro steps, 
depends on the ratio of the widths of the and n re- 
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FIG. 3: (a) Current-voltage characteristics showing rf- 
induced Shapiro steps both at the usual voltages nhf r f/2e 
and at half-integer values nhf r f /4e. (b) Temperature de- 
pendence of the maximum (power-optimized) critical current 
steps, showing that the integer steps scale with the junction 
critical current whereas the half-integer steps only occur at 
temperatures near the minimum in the critical current. 
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FIG. 4: Phase diagram mapping the region of wq/Xjo and 
Wtt/Xjtt values for which a spontaneous circulating current 
flows around the 0-7r step edge. In the shaded region, the 
phase drop across the junction <j> and current density J vary 
spatially along the junction width. Outside this region, the 
phase is uniform with value or n and no spontaneous cur- 
rents flow in the junction. 



gions of the junction wq and u) T to their correspond- 
ing Josephson penetration depths Ajo and Aj w 
Here, Xj = (h/2e/j,odmJ c ) 1 ^ 2 , where the magnetic bar- 
rier thickness d m = 2X + d depends on the penetration 
depth of the superconductor A and the barrier thickness 
d. The phase diagram denoting the regimes of uniform 
phase and spontaneous current flow is shown in Fig. 0] 
Regimes with uniform phase differences <fi — or <f) — n 
are separated by a narrow region in which the phase dif- 
ferences vary across the junction due to spontaneous cir- 
culating currents. The circulating currents onset along 
lines defined by Xj n tanh^o/Ajo) = A jn tarj/Wn- IXj-A 
and Ajo tanh(w7r/A,/ w ) = Xj^ tan(wo/Ajo) |XJJ| . Also 
shown is the path in the phase diagram followed by the 
non-uniform junction described in Fig. [^a), for which 
d m = 100 nm, as the temperature is varied from K to 
2.1 K, the temperature at which the wide segment of the 
junction crosses into the state. We see in the inset that 
spontaneous currents occur only in a narrow tempera- 
ture range from 1.54-1.56 K. This is comparable to the 
temperature range from 1.50-1.56 K in which we observe 
half-integer Shapiro steps, suggesting that they have the 
same origin. Other non-uniform junctions that we have 
measured either had smaller areas so that wo and 
were smaller, or smaller barrier thickness steps so that 
Ajo an d Aj w were larger. In either case, the tempera- 
ture range of spontaneous currents was predicted to be 
immeasurably small (< 1 mK), explaining why we were 
not able to observe half-integer Shapiro steps or a finite 
minimum critical current in these junctions. 

Previously, we directly measured the current-phase re- 
lation of uniform SFS junctions with thicker barriers d « 
22 nm, close to the second node of the I c (d) dependence, 



and did not observe any signature of a sin(2</>) component 
in the Josephson tunneling Q • This is in agreement with 
microscopic predictions that the second-order Josephson 
tunneling probability falls off more rapidly with increas- 
ing barrier thickness than the first-order Josephson effect 
Thus, to determine definitively whether or 
not a sin(20) component is present in the current-phase 
relation, uniform SFS Josephson junctions with barrier 
thicknesses d ks 11 nm close to the first node of I c (d) 
should be studied. Simultaneous observation of Shapiro 
steps at half-integer Josephson voltages and period dou- 
bling in the current-phase relation in uniform SFS junc- 
tions is required to reliably verify the presence of a mi- 
croscopic sin(2</>) component. 
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